Narrow resonances and black-hole-like absorption in a non-black-hole metric 
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A massive body with the Schwarzschild interior metric (perfect fluid of constant density) develops 
a pressure singularity at the origin when the radius of the body R approaches 9r s /8, where r 3 is the 
Schwarzschild radius. We show that a quantum scalar particle scattered in this gravitational field 
possesses a dense spectrum of narrow resonances. Their density and lifetimes tend to infinity in the 
limit R — s> 9r s /8, and we determine the cross section of the particle capture into these quasibound 
states. Therefore, a body that is not a black hole demonstrates black-hole-like absorption. 

PACS numbers: 04.62.+V, 04.70.Dy, 04.70.-s 



I. INTRODUCTION 

We found in our previous work [l[ that the scattering 
of a massless scalar particle in a near-black-hole met- 
ric of a massive body whose radius R slightly exceeds 
the Schwarzschild radius r s , is characterized by a dense 
spectrum of resonances. For R — > r s both the resonance 
energy spacing D and their width 7 tend to zero, while 
their ratio remains finite (7/-D ~ 2e 2 r 2 /iT) and tends to 
zero for small energies e. (We use units where % = c = 1.) 
This allowed us to define the cross section for particle 
capture into these long-lived states in the spirit of the 
optical model [2[ , by averaging over a small energy inter- 
val containing many resonances. Note that this capture 
emerges in a purely potential scattering problem, without 
any absorption introduced a priori. Somewhat unexpect- 
edly, the capture cross section turned out to be equal to 
the cross section obtained by assuming total absorption 
at the event horizon. 

In particular, in the zero-energy limit our result coin- 
cides with Unruh's absorption cross section a a = 4irr 2 
for a black hole 043- This shows that a non-singular 
static metric can acquire black-hole properties prior to 
the actual formation of the black hole. 

It is interesting to see whether this resonance absorp- 
tion phenomenon is specific to near-black-hole metrics, 
or if there are other instances of similar behaviour. For 
example, consider a massive star modelled as an incom- 
pressible fluid sphere. The interior of such a body is 
described by the Schwarzschild interior metric 10]. This 
metric develops a pressure singularity at the origin for 
R -> 9r s /8, far form the black-hole limit. In Ref. [ll[ 
we examined the spectrum of discrete bound states for a 
massive scalar particle in the gravitational field described 
by this metric. It was found that near the singularity 
this spectrum acquires some black-hole-like features. In 
particular, all the levels with finite principal quantum 
number n collapse towards zero energy (binding energy 
is —mc 2 ), i.e., the spectral density becomes infinite, sim- 



This gives us a motivation to search for narrow reso- 
nances in the Schwarzschild interior metric. Indeed, in 
the present work we find such resonances in the limit 
R — > 9r s /8. However, in contrast to the black-hole-like 
case considered in Ref. |l|, the effective potential pro- 
duced by the Schwarzschild interior metric possesses a 
barrier near the boundary, and the long resonance life- 
times are due to the particle having to tunnel through 
this barrier to the surface of the body. 

Similar to the black- hole- like behaviour for R — > r s , the 
resonance width 7 and energy spacing D tend to zero for 
the Schwarzschild interior metric in the limit R — > 9r s /8, 
while their ratio j/D remains finite and tends to zero 
at small energies e. One can thus define the "optical" 
capture cross section by averaging over a small energy 
interval containing many resonances. In contrast to the 
near-black-hole case, where the cross section remains fi- 
nite, the absorption cross section for the Schwarzschild 
interior metric tends to zero at zero energy. However, 
this metric still possesses black-hole-like absorption for 
nonzero energies. This means that such an object may 
gravitationally absorb particles in a black-hole-like man- 
ner despite the fact that its metric never approaches that 
of a black hole. 

In what follows we analyse the massless scalar particle 
scattering problem for the Schwarzschild interior metric 
both numerically and analytically. Although the latter 
treatment of the problem is only approximate, it con- 
firms all the important features of the emerging resonant 
scattering picture. 



II. RADIAL KLEIN-GORDON EQUATION 

The Klein- Gordon equation for a scalar particle of 
mass m in a curved space-time with the metric g^ is 



-gm 2 ^ 



0. 



(1) 



ilar to that of a near-black- hole metric (cf . Ref. [12| ) . 



For a particle of energy e in a spherically symmetric 
field we seek solution of Eq. (JXJ) in the form ^>(x) = 



2 



e~ let ip(r)Yi m (9, <p). Considering for simplicity the case of 
a massless particle with zero angular momentum (Z = 0) 
and the metric of the form 



ds 2 = a{r)dt 2 - b{r)dr 2 - r 2 dtt 2 , 



(2) 



where dfl 2 = d6 2 + sin 2 6dip 2 , the radial wave equation is 
given by 
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where h(r) = y/a(r)/b(r). By transforming the radial 
function as ip(r) = r _1 <j){r) / ' \J 'h(r) , Eq. can be cast 
in the following Schrodinger-like form, 



(4) 



This equation is convenient for deriving its semiclas- 
sical (WKB) solution. For small h(r) (near r = 0, see 
below) the first term in braces dominates and this solu- 
tion is 
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In fact, it is easy to verify that this solution is exact if 
the —h'/rh term in braces in Eq. ((4]) is neglected. 



III. PHASESHIFT FOR THE SCHWARZSCHILD 
INTERIOR METRIC 

Outside a spherically symmetric, non-rotating body 
of mass M and radius R the metric is given by the 
Schwarzschild solution 



ds 2 = { 1 - — )dt 2 - 1 - 

r / V r 



dr 2 ~r 2 dfl 2 , (6) 



where r s = 2GM is the Schwarzschild radius of the body 
and G is the gravitational constant. Hence, for r > R 
Eq. (J3J takes the form 
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r^(r)=0. (7) 



The Schwarzschild interior solution describes a static 
spherical mass of incompressible perfect fluid with con- 
stant density, and is given by Eq. ([2]) with 



6(r) = 1 - 
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so that 
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FIG. 1: The short-range phaseshift S obtained numerically 
as a function of e = eR (solid line) and the semiclassical 
phase, Eq. (|13p , at the classical turning point (thick dashed 
line) for £ = 8/9 - r s /R = 0.00014. The vertical dashed line 
corresponds to the sixth resonance at e w 0.1080. 



This metric is valid for r s < 8R/9 but develops a singu- 
larity as a(0) vanishes for r s = 8R/9 13]. 

Using h(r) from Eq. ((5J) in Eq. gives the radial wave 
equation for r < R. For r s < 8R/9 this metric is smooth 
at the origin. Hence, the regular solution for Z = must 
satisfy the boundary conditions i/j'(0) = 0, ^(0) 7^ 0. 
The value of ip(0) only affects the normalization of the 
wavefunction and we set tp(0) = 1. We solve the interior 
equation numerically using Mathematica [i~4| . 

This solution supplies the boundary condition for the 
exterior wavefunction at the surface of the body r — R. 
(R — 1 is used in the numerical calculations.) Equation 
(O is then integrated outwards to large distances r ^> r s . 
In this asymptotic region Eq. ([7]) takes the form of the 
nonrelativistic Shrodinger equation for a particle with 
momentum e and unit mass in the attractive Coulomb 
potential with charge Z = —r s e 2 . Here the wavefunction 
can be matched with the asymptotic Coulomb solution 
2], 



i/j(r) oc sin[er — (Z/e) \n2er + 5c + 5] 



(9) 



where 5c = argT(l + iZ/e) is the Coulomb phaseshift, 
to find the short-range phaseshift 5. The latter is de- 
termined almost exclusively by the interior metric, and 
carries important information about the behaviour of the 
wavefunction at r < R. 

Unlike 5c which is small and has a weak dependence 
on the energy, the phaseshift 5 becomes large when r s 
approaches 8R/9, i.e., for 



£ = 8/9 - r s /R < 1. 



(10) 



This phaseshift also has a strong dependence on the en- 
ergy of the particle, as shown by the solid line in Fig. [1] 
Similar to the case of a near-black-hole metric consid- 
ered in Ref. [l[, this phaseshift goes through many steps 
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FIG. 2: The radial wavefunction at the energy e = 0.1080 of 
the sixth resonance for £ = 0.00014 (solid curve). The dashed 
curve is the analytical approximation for the wavefunction 
given by Eqn. ()16[) for the energy e„ (n = 6) from Eqn. (|19|) . 
The vertical lines indicate the two classical turning points 
defined by p(r) = 0. Regions I and III are classically allowed, 
while II is classically forbidden. 



In the case of the Schwarzschild interior metric, h(r) is 
given by Eq. ©. For £ <C 1 and r -C R (in practice r < 
0.5R is sufficient) h(r) can be approximated as follows, 
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(12) 



This shows that for small r, where the first term in braces 
in Eq. (fTT|) dominates, the particle has a large classi- 
cal momentum (e.g., p(0) ~ 4e/9£). This corresponds 
to the motion in a deep classically allowed region near 
the origin. Analysis of Eqs. ((8]) and (fTT|) shows that for 
£ <C 1 there is a second classically allowed region near the 
boundary with a broad, nonsemiclassical potential bar- 
rier in between, cf. Fig. [2] These regions are separated 
by two classical turning points r\ and given by the 
roots of p(r) = 0. 

The above picture explains the origins of the dense 
spectrum of narrow resonances observed in the limit 
£ — » 0. To estimate the resonance widths analytically 
we require an expression for the scattering matrix, which 
is constructed using the procedure outlined in the follow- 
ing subsections. 



of the size ir. Each of these steps corresponds to a reso- 
nance, i.e., a long-lived quasibound states of the projec- 
tile and the target. For energies corresponding to the 
midpoints of the steps (where the derivative dS/de is 
largest) the magnitude of the wavefunction ip(r) inside 
the body (r < R) is much greater than outside. This is 
a signature of a quasibound state. The wavefunction at 
one of the resonances is shown in Fig. [2] 

The resonances observed in the phaseshift S are char- 
terised by their energies e n and widths j n , which can 
be determined numerically by fitting the "steps" in 6 
by arctan[2(e — e n )/7n] + const. @], or analytically (see 
Sec. ITV)) . We show below (see Figs. [3] and @| that 
as the Schwarzschild radius r s tends to 8R/9, i.e., for 
£ — » 0, both the energies and the widths of the reso- 
nances tend to zero (i.e., their lifetimes tend to infinity). 
This means that in this limit the massive body with the 
Schwarzschild interior metric develops absorption prop- 
erties usually typical of black holes. 



1. Region I: origin to barrier 

In this region (0 < r -C r\) the e 2 /h 2 term in Eq. (fill 
dominates near the origin for £ -C 1. Dropping the last 
term in Eq. (fT2|) we obtain the semiclassical phase of the 
wavefunction <f>(r), Eq. (|S"|), 
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where e = eR and p — r/R are the scaled energy and 
radial coordinate, respectively. The wavefunction in this 
region is then 
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where C depends on normalization. 



IV. ENERGIES AND WIDTHS OF THE 
RESONANCES 



2. Region II: suppression by barrier 



A. Schwarzschild interior 

When the wave equation is written in the Schrodinger- 
like form (UJ), the coefficient in braces playes the role of 
an effective potential for the motion a particle with the 
classical momentum 
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In the intermediate range, < p < 0.5, the square 
of the classical momentum ([TT]) simplifies to 



P 2 (r) 
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4r 4 



(15) 



In this case Eq. (|4]) has an exact analytical solution, 
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FIG. 3: Energies of the n = 4 resonance given by Eq. (|19p 
(solid line) and obtained numerically by fitting the resonance 
jumps of the phaseshift S (solid circles) as a function of r a /R. 



where $ and A are constants determined by matching to 
the solution (IT4l) to the left of the first classical turning 
point p x = 9e/\/8 [from Eq. (fPoj) ]. 

For < P < Pi both Eq. (JUJ) and Eq. ([T6j are valid 
(the former is valid at the origin while the latter is not). 
Matching is done by first expanding the arctan term in 
Eq. (|14[) for p/v 7 ? 3> lj and neglecting the arctan term 
in Eq. (fTo| . which is justified for 2p/9e <C 1. This yields 
the phase 



$ = 



2£ 



(17) 



When the energy is on resonance, the wavefunction 
decreases under the barrier (Region II in Fig. [2]). Con- 
sidering Eq. (|16l) for p 3> px (but keeping p < 0.5, which 
is possible for low energies e <C 1), and expanding the 
arctan function for large arguments, we have 



bnifi) ~ -2Ap 2 sin 
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This function represents a decreasing solution <j)\i{p) oc 
only if we require Tre/i/2^ — 7r/2 — nir (n = 1,2,...), 
which gives the resonance energies 



(19) 



2£(n + l/2). 



Hence the resonances form an equispaced spectrum sim- 
ilar to that of a harmonic oscillator (starting from n = 1 
though). This is in agreement with the behaviour of the 
phaseshift in Fig. [1] 

According to Eq. (fl"9|). for £ — > the energies of all 
resonances tend to zero. This behaviour is shown in Fig. [3] 
which compares the analytical expression for e„ with the 
values obtained numerically from the phaseshift <5. 



3. Region III: barrier to boundary 

To calculate the S-matrix and determine the resonance 
widths analytically we need to match the interior and 
exterior solutions at the boundary r — R. Formally, the 
effective potential in Eq. Q contains a second classical 
turning point r2, which defines the outer classically al- 
lowed region, r^ < r < R. In this region the classical 
momentum remains small and the wavefunction changes 
little between r^ and R (see Fig. [5]). Indeed, it is easy 
to check that the classical momentum takes its largest 
values at the boundary r = R, and for £ —> one 
has p(R) = V81e 2 + 58i? -2 . This momentum is much 
smaller than the momenta at the origin (see Sec. IIV A ip 
or at small r [see Eq. (|15[)]. 

Therefore in generating an approximate analytical so- 
lution over the entire interior region we assume that the 
major contributions arise from regions I and II. Changes 
in the wavefunction due to the potential in region III give 
rise to higher-order effects which we calculate as correc- 
tions at the end of this section. Thus we formally extend 
the solution in Eq. (1161) to the boundary p — 1 where 
it is matched to an appropriate solution of the exterior 
equation as discussed below. 



B. Schwarzchild exterior solution 

For low-energy scattering, eR <C 1, the last term in the 
exterior equation ([7]) can be neglected near the boundary. 
Hence, for r ~ i?, the solution of Eq. is 



ip(r) — a In 
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where a and /3 are constants determined by joining this 
solution with the interior solution at r = R (see below). 
Following the matching procedure outlined in Ref. [15[ 
the low-energy s-wave scattering matrix is found as So = 
S e 2iSc , where 
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(PI a) - isr s C 2 



(21) 



is the short-range part of the scattering matrix, and C 2 = 
2irer s /[l ~ exp(— 27rer s )] is a function that characterizes 
the long-range Coulomb tail of the potential • 

In general, resonances can be found as poles of the S'- 
matrix in the complex energy plane at e = e n — i~f n /2 0. 
On the real energy axis they correspond to the points 
where the short-range phaseshift S passes through 7r(n — 
1/2). At these energies So = — 1, which requires (3/a = 0. 
This means that near the resonances one can use 



p/a~{p/a)'{e 



(22) 



where the prime denotes the derivative with respect to e 
taken at the resonance. Using this expression and looking 
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for the poles of So, Eq. (j2"Tj) . gives the resonance width 
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To find the ratio (3/ a and (/?/«)' we match the loga- 
rithmic derivatives of the interior and exterior wavefunc- 
tions at the boundary. Using the relation between the 
solutions ip( r ) ancl 4>( r ) (Sec. [TTJ) for the interior met- 
ric, Eq. ([8]), the logarithmic derivative (evaluated at the 
boundary) is 
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where we used Eq. ([16]) to obtain the last line. 

The exterior logarithmic derivative at the boundary for 
r s -» 8/9-R is found from Eq. (f!HI)) as 
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FIG. 4: Width of the n — 4 resonance found numerically 
by fitting the phaseshift 8 (solid circles), and estimated from 
Eqs. (Tl9|) and ([28} (solid line), as a function of r s /R. The 
inset shows the same data as a function of £ on a double 
logarithmic scale. 



ABSORPTION CROSS SECTION 
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Setting the logarithmic derivatives from Eq. (j!M]) and 
Eq. (p5)l equal, gives 
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When finding (/3/a)' to estimate 7„ from Eqn. (f23| . we 
can use the fact that at (or very near to) the resonance, 
(3/a = 0. This allows one to simplify the answer and 
obtain 



1; 



(27) 



with the constant B given approximately by B = 2.7 '4. 

Thus we see that similar to the resonance energies e„, 
Eq. (fT9|) . the resonance widths also decrease with £, i.e., 
as the metric singularity is appoached. Taking Eq. p^|) 
into account we see from Eq. (p?T|) that for a fixed reso- 
nance number n, j n cx £ 5 / 2 . This dependence is shown 
on the inset of Fig. [4] 

The above derivation of the width underestimates the 
width (through the value of B) by a factor r\jR 2 ~ 0.1, 
since the correct wavefunction does not decrease through 
region III as quickly as that given by Eq. (fT5|) (see Fig. [5]). 
Restoring this factor and evaluating B numerically gives 



In the limit r s — > 8R/9 the short-range phaseshift 
and the energy density of the resonances tend to infinity. 
Hence, neither d, nor the resonance energies and widths, 
e n are 7„ (which tends to zero) retain much physical 
meaning. However, as in the near-black-hole metric case 
[l|, the cross section for capture into these resonances, 
i.e., the effective absorption cross section, is well defined. 

Equation (fl9|) shows the level spacing between the res- 
onances is given by 



D = e n+1 - e„ = V^/R- 



(29) 



Assuming a finite energy resolution As ^> D and given 
that 7„ <C D, one can introduce the cross section for 
capture of the particle into these long-lived states. Ac- 
cording to the optical- model considerations Q, this ab- 
sorption cross section is given by 
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Substituting the expression for j n , 
above expression yields 
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Unlike the low-energy absorption cross section for the 
black holes, a a = iirr 2 , the cross section (|3"Tj) vanishes 
at zero energy. However, for non-zero incident energies 
the capture capture cross section for massless particles is 
finite. 



In 



28^e 4 n R 3 , 



(28) 



VI. CONCLUSIONS 



in good agreement with the numerical calculation, as 
shown by Fig. 



We have considered the problem of scattering of mass- 
less scalar particles from a spherically symmetric con- 
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stant density fluid sphere described by the Schwarzschild 
interior metric. We find that despite the black hole limit 
being unattainable for such metric, a dense spectrum of 
narrow resonances emerges in the limit r s — > 8R/9 in 
which a singularity develops in the metric at the origin. 



This phenomenon gives rise to a nonzero capture cross 
section of massless particles for nonzero energis. This 
implies that an object that is not a black hole may grav- 
itationally absorb particles making it appear black hole- 
like, albeit with a different capture cross section. 
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